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1 Introduction 



This article deals with a stochastic Poisson equation on a bounded domain DcK', 
with arbitrary dimension k > 1, driven by a fractional Brownian field B H , with 
H = (Hi, . . . ,H k ) G [\, l[ k . We prove a theorem on existence and uniqueness of 
solution, we study the properties of its sample paths and finally, we give a numerical 
scheme based on lattice approximations and prove its convergence on a functional 
space with some explicit rate. 

The equation is described as follows: 



(1) 



Au(x) - f(u(x)) = g(x) + B H (x), x G D, 
u(x) = 0, x G 3D. 

We assume that / has a decomposition f = fi + f 2 , with / l5 f 2 : ffi — > R satisfying 

(fl) /i is continuous, non-decreasing and sup,j. 6K < M, 

(f2) /2 is Lipschitz with small Lipschitz constant L. 

The function g : D — > R is measurable and satisfies some integrability conditions. 
The stochastic character of the equation comes from B H (x), which denotes the 
formal derivative of a fractional Brownian field. 

We give a rigorous meaning to (1) by means of a mild formulation, as it is pretty 
usual in the SPDEs literature. For this, we recall that if k > 2, the Green function of 
the deterministic Poisson equation on a bounded domain D with smooth boundary 
is given by 

G k D (x,y) = G k (x,y) -E x (G k (B T ,y)), x,y G D, (2) 



with 



G k (x,y) = C k < 



\og\x-y\, k — 2, 
\x-y\ 2 - k , k>3. 



Here C 2 = C k = k ( 2 -k)uj k ^ OT k > 3, where cu k denotes the volume of the unit ball 
in M, k , and B T is the random variable obtained by stopping a A;-dimensional Brownian 
motion starting at x at its first exit time of D (see for instance [9] and also [7]). 
For k — 1, G l (x,y) = C\\x — y\ (see for instance [17], pg. 16). The expression for 
G]j(x, y) depends of the domain D. If D =]0, 1[ then G l D (x, y) = (xAy) — xy, where 
"A" denotes the infimum (see [7] pg. 258). 

By a solution to (1) we mean a stochastic process u = {u(x), x G D} satisfying 

u(x)= f G k D (x,y)f(u(y))dy+ f G k D (x,y)g(y)dy + / G k D (x,y)dB H (y). (3) 
Jd Jd Jd 

For a similar SPDE in dimensions k = 1, 2, 3, driven by a standard Wiener field W, 
different problems have been studied so far. For instance, existence and uniqueness 
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of solution has been proved in [4] using the classical theory of equations defined by 
monotone operators (see [27]); the Markov field property has been investigated in [5] 
and [6], and numerical approximations have been given in [11]. Let us remark that 
for k < 3 the stochastic convolution J D G k D (x,y)dW(y) is well defined as a Wiener 
integral, because G k D (x, ■) is square integrable. For k > 4, this property is not true 
anymore. 

For k > 4, a SPDE of the same type than (3), driven by a Gaussian stationary 
process F with an absolutely continuous covariance measure, but possibly having 
singularities, has been studied in [20], extending the results of [11] to higher dimen- 
sions. The authors combine conditions on deterministic functions and covariance 
densities derived from Young's type inequalities and provide a definition of an in- 
tegral with respect to the random field F, and thereby a suitable meaning of the 
stochastic convolution J D G k D (x, y)dF(y). As regards the approximation scheme, the 
approach for k < 3 based on a Fourier series expansion of G k D (x,y) is not appro- 
priate. In fact, as has been already mentioned, G k D (x, •) is not square integrable 
and therefore this function and its Fourier series may not coincide. Instead, a more 
sophisticated procedure involving a smoothing of the Green function combined with 
its Fourier series expansion has been considered. 

The main reason for considering correlated noises in high dimensions is to com- 
pensate the irregularity of the Green function. However, it may also be a rea- 
sonable choice when modeling phenomena where the stochastic input shows some 
dependence. 

With the increasing attention devoted to fractional Brownian motion in recent 
years, the study of different type of problems on SPDEs driven by fractional noise is 
being more present in the mathematical literature. We refer the reader to [26] for an 
extensive list of references on the subject, including some motivating aspects from 
other disciplines. Generically, these SPDEs are appropriate to model phenomena 
showing up either persistence (for example, in hydrology or finance) or intermittency 
(as turbulence). At the best of our knowledge, developments on this topic refer so 
far mainly to parabolic and hyperbolic SPDEs, the elliptic case being less explored. 

A particular version of Equation (3) with null functions / and g appears in [14] 
(for Hi e]|, l[ k ) and in [15] (for Hi e]0, l[ fc ). In both references, the authors apply 
white noise analysis to give a meaning to the solution u(x) = j D G k D (x,y)dB H (y) 
in the sense of distributions. Conditions on Hi ensuring the existence of an L 2 (D)- 
valued solution are given. In comparison with these references, our analysis of (3) 
allows a monotone nonlinearity f(u) and a free term given by g, as in [11] and [20]. 

Very recently, in [19], a new and very promising approach to stochastic Poisson 
equations based on Wiener chaos solutions on weighted spaces has been given. Using 
Malliavin calculus and a formulation of the stochastic convolution by means of the 
Wick product, the authors prove existence and uniqueness of solution for elliptic 
equations allowing some type of nonlinearities in the noisy term, modelled by a 
general spatial Gaussian process. Furthermore, in [16] numerical approximations 
based on a finite element procedure are provided. 

The content of the paper is as follows. In Section 2, we give some preliminaries 
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on the fractional Brownian field B H , when H = (Hi, . . . , Hk) G [|, l[ fe . We combine 
ideas from [14] and [15] (see also [8]) with some results from [23] and [22] to give a 
moving average representation of B H in terms of an standard Wiener field. Then, 
we identify a suitable L p -space with mixed norm of deterministic functions which 
can be integrated against B H . These spaces are related with the reproducing kernel 
Hilbert space of B H by means of Hardy-Sobolev's inequality. Section 3 is devoted to 
the study of the stochastic convolution of the Poisson kernel (2). We give a sufficient 
condition on the Hurst parameter H ensuring the integrability of the Poisson kernel 
with respect to B H , according to the result proved in Section 1. We also give some 
probabilistic properties of the stochastic convolution and prove the Holder continuity 
of its sample paths (Theorem 3.3). These ingredients, though of own interest, are 
meant to provide a rigorous meaning to Equation (3) in any dimension k > 1. In 
Section 4, we give a theorem on existence and uniqueness of solution of Equation (3) 
on the space of continuous functions vanishing at the boundary dD; we also prove 
Holder continuity of the sample paths of the solution. Finally, Section 5 is devoted 
to numerical approximations of (3). We consider the domain D =]0, l[ k and use the 
approach of [11] for k < 3 and that of [20] when k > 4. With an appropriate choice 
of the functional spaces, we give the rate of convergence. For k < 3 we find the 
same rate as for the Brownian case, while in dimensions k > 4, the rate depends on 
the regularity of the noise, as may be expected. 

Throughout the paper, we shall denote by ch any positive constant depending 
on the Hurst parameter H = (Hi, . . .Hk), k > 1, independently of its particular 
value, and by C any positive, finite constant. 

2 Preliminaries 

Let H = (Hi, ...,H k ) G]0, l[ fc . A fractional Brownian field (ffif) on R k with Hurst 
parameter H is a Gaussian stochastic processes B H = {B H (x), x G R k } with zero 
mean and covariance function given by 

k 

R H (x,y)=E(B H (x)B H (y)) = ]J R Ht (x l ,y t ), (4) 

i=i 

where 

Such a process has been introduced and considered in relation with different prob- 
lems in [12], [13], [14], [15]. 

In this article, we consider values of the Hurst parameter H = (Hi, . . . , H k ) G 
[|, l[ fe . Our first goal is to define a stochastic convolution for the Poisson kernel with 
respect to B H . For this, we shall identify a suitable class of deterministic functions 
/ : R k -> R for which 

JR k 
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is a well denned random variable. As in the one parameter case, it will be useful to 
have a moving average type representation of the process B H in terms of a standard 
Brownian field on R fe . We shall prove such a representation owing ideas from [15] 
but considering the framework of [23] and [22] (see also Lemma 1.20.10 in [21]). 

We start by introducing some notation. On R fc we consider the usual partial 
order defined coordinatewise and denote by x = (x±, . . . ,X)~) a generic element in 
this space. For x,y G R fc satisfying x < y, we set t[ x>y ](rj) = Y\ k =1 \ XuVi \(r]i). We 
shall denote by £ the set of elementary functions on R fe , that is functions of the 
form 

lo 

<p(v) = 

i=i 

with <pi G R and disjoint rectangles [x l , y l ], I — 1, . . . , lo- For <p G £, we define 



l( V ) = J2^B H ([x l ,y 1 ]), 



i=i 



where B H ([x l ,y l ]j denotes the increment of B H on the rectangle [x l ,y l ] in the 
sense of A;-dimensional distribution functions. For any x G R fc , we set l[o tX ](v) = 
nil %,xi](»7i), where 

[ 1, Vi e [0,2*], 
%,xi](»7i) = S -1, ffc e [xi.O], 
[ 0, otherwise. 

Then on £ , we introduce the inner product (•, -)^h derived from the covariance 
structure of B H given by 

(%),*], %>,»])«» = Rh{x,v). (5) 

The closure of £ with respect to the norm || • \\ h h will be denoted by 7i H . 

For each i — 1, . . . , k, we define the linear operator acting on functions y? G £ as 
follows: 



(v) 
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CH i J R v{rii,---,Vi-i,Ui,r] i+1 ,...,ri k )(u i -?&)+' 2 c?Mi, #;G]§,1[, 



where are constants depending only on (see (2.9) in [15] for its explicit value). 
For its further use, we introduce the sets = {i = l,...,k : Hi > |}, 

= {i = 1, . . . , k : Hi = ^} , and denote by c^, c^ = \ the cardinals of and 
C^, respectively. We notice that 

( k h^) (n) = f • • • ■ ■ ■ ,Vk)vU du i)> ( 6 ) 



where 



CH 4 (iii - ??i)+ l 2 *^, z G C (>) , 
S Vi (dui), ieC(=\ 



and 8 Vi denotes the Dirac measure at rji. 
By iteration, for k > 2 we define 



Notice that for a function <p = ®f =1 </?i, 

(^ (fc V) (v) = n ( K k<fi) (vd- (8) 
i=i 

Consider the one-dimensional case (k — 1). The Mandelbrot- van Ness represen- 
tation establishes that for any H g]0, 1[, 



B H (x) = / K H (x,y)W(dy), 



where 



K H (x, y) = c H {(x - y)J _ * " (-tf)?"*} - ( 9 ) 
and is a standard Brownian motion on the real line (see [25], pg. 320). Thus, 

K H (x, v)K H (y, v)dv = R H (x, y). 
A simple computation yields [k* h 1[o, x ]) (y) — K H (x,y), for H > \. Hence, 
T (%,,„]) := = f (K* H 1 [0 , X] ) (y)W(dy). 

With the definition (6) this representation also holds for H — \. Thus, it is clear that 
the mapping l [0jX] ^ J R K* H (t [0 ^ (y)W(dy) is an isometry between (£,(•, -) h h) 
and L 2 (Q). That means, the operator K* H is an isometry between (S, (•, and 
L 2 (R). Hence, K* H can be extended to the Hilbert space 7i H . Otherwise stated, 
we can extend the definition of X from £ to 7i H , and therefore define Z(y?) = 
/rC-^h^X?/)^^?/) as an £ 2 (^)-valued random variable. 

In the multidimensional case we prove similar results. We first recall a definition 
and introduce some notation. A stochastic process {W(x),x G R fe } is termed a 
standard Wiener field on R fc if it is Gaussian, with mean zero and covariance given 
by E (W(x)W(y)) = x Ay, where x A y := Ui=i( x i A y*), and 



XiAyi, x i} yi>0, 
(-Xi)A(-yi), Xi,yi<0, 
0, otherwise. 



For a function ip : R — > R, we define </?(iti, . . . , u^, V\, . . . , t^) = ip(wi, . . . , it^), with 

the set of 



u>i = Ui, if i G ) and Wj = ^ if i G C^. Then we denote by 7i 
functions </? : R fc — > R such that 

|2ff;-2 

yuuiUViiii^iii — L)\Ui — Vi 

iec(=) iec(» 

X ^(uj, . . . ,M fc )||(^(u 1 , . . . ,U fc ; Vx, . . .,v fc )| < +oo 



H 



II / r2c(>) II {diHdViH&Hi - l)\ Ui 
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Proposition 2.1 Set Kh{x, y) = Y\ k =1 Kn^Xi, yi), with Kh^ i — 1, • • • , k, defined 
in (9). Then: 

1. There exists a standard Wiener field W on R fe such that 

B H (x) = j Rk K H (x,y)W(dy) = (K*/ k % x] ) (y)W(dy). (10) 

2. For any <p G H H 

Xitp) = J Rk (y)dW(y) (11) 

defines a random variable in L 2 (Q). 

3. For any ipi,<fi2 £ 'H H the following isometry formula holds: 

E(lMl(<p 2 )) = (^ (fe Vi) (V) {K*/ k ^ 2 ) (y)dy 

= L=, n ^L» n {duMm - i)\ Ut - v t \^- 2 ) 

m ied=) JR ied» 

x (p 1 (ui,...,u k )ip 2 (ui,---,u k ;vi,...,v k ). (12) 
In particular, if Hi e]|, 1[ /or any z = 1, . . . , fc, 

A; 

E(J(^)J(^ 2 ))= / ^( M )^ 2 (^f[(^(2^-l)|^-^| 2 ^- 2 )^. (13) 



i=l 



Proo/ To prove the existence of W, we follow the arguments of [22], pg. 279, for 
k — 1, which extend easily to any k > 1, as follows. 

For each z = 1, . . . , k, such that ifj > |, the action of the kernel ATjy. on elemen- 
tary functions can be expressed in terms of a fractional integral. More precisely, if 

<p G £, K* Hi tp = c H J H *~*(p, where jf* = r(g--i) - x) H ^^(p(y)dy (see 

section 3.2 in [23]). Then, for ^ in the image of K* H and by considering the frac- 

H — - l 

tional derivative defined by D_ l 2 if)(x) = r ,t]_„ , (ip(x) — ij){x + y)) y~?~ Hi dy 

H ■ — - 

we define = c ~kiD- 2 ifj, which by the rules of fractional calculus is seen to be 

the inverse operator of Kjj . For — |, Q* H . is defined to be the identity operator. 

Set W(x) = B H (n*U Q* H J[ , Xl ])- The process W = {w(x), x G R k } is a stan- 
dard Wiener field on M fc , and the stochastic field B H has the integral representation 

B H (x) = f K H (x,y)dW(y). (14) 

Indeed, set Q^ (fc) l [0tX] = U-=i W Usin S ( 8 )' we have ( K H ik) ° Q*^) h,*] = 
l[ 0jX ]. Thus, for any x,y E M k , 

E(W(x)W(y)) = {Q* H {k H %xh Q* H ik) \z, y] ) HH 
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By construction, the operator K*^^ is an isometry from S into L 2 (IR fc ) that can be 
extended to the Hilbert space H H . Therefore, one can define T(ip) for any (p e H H 
by means of (11). 

We now prove (12). By the very definition of K*^ (see (6)) and by applying 
Fubini's theorem, we obtain 



k (y) K (fc V 2 ) (y)dy 

L^iL^ n n ^(*o 



X 



iec(=) iec(» 

^(v) n n $( dv i) 

iec(=) iec(» 



»ec(=) 



L(» n du * dv 4 L» n ^4>*-y*)f 



X 

x y>i(ui, . . .,u fc )^ 2 (Mi, . . . . . .,«*;). 

From this and the identity 

4 / dy t { Ul - yi )+~'k^ ~ yi)+~~ 2 = H t (2H t - 1)K - v t \ 2H >- 2 

(see [10], page 404), (12) follows. 

Finally, if C» = {1, . . . , k}, (12) reads (13). This ends the proof of the Propo- 
sition. 

■ 

It is well known that for real functions tp, ip and H e]|, 1[, 

JJ R Hv)\\m\\v - e\ 2H - 2 dvde < MMU (M) IHU (R) , (is) 

with a positive constant bn- Indeed, this follows from Hardy-Littlewood-Sobolev's 
inequality (see inequality (1), page 321, in [3]). By considering functions <p\,<f2 '■ 
18L k — > IR and applying recursively this inequality for indices i E and Schwarz 

inequality for i e C^ = \ we obtain 



l)\Ui-V, 

iec(=) iec(» 

X ]</?i(mi, . . . , Ufc)||^2(«l, • • • , Uk, Vi, . . . , V k 



\2Hi-2^ 



< C//i HtllVlll —L- -J_ \\f2W _L- _L- , (16) 



where for pi G [1, 00], i = 1, . . . , k, 

— Vk 

|LFi.-.Pfc(ffifc) = f / (/•••(/ V)k)\ VX dT)\ J P1 ••• j"*" 1 ^* 
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Let us denote by L Pl '-> Pk (R k ) the space of measurable functions h defined on R 
with ||/i||i,pi,-,Pfc(Rfc) < oo. Such spaces are termed LP spaces with mixed norm. For 
details we refer the reader to [3] and also [1]. In particular, if Hi e]|, 1[ for any 
i — 1, . . . , k, a proof of (16) is given in page 322 of [3], but it is easy to extend the 
result allowing the value ifj — | for some indices i. The preceding discussion yields 



i i 



L H i'-' H k(R k ) c n H c n H . (17) 

For its further use, we remark that for any p > sup ie ^ k yPi, and every mea- 
surable function h with bounded support O contained in R fc , 

\\h\\ LPl ,-, PHO ) < C\\h\\ LP{0 ), (18) 

with a constant C depending only on O. Indeed, this follows by applying recursively 
Holder's inequality with pp^ 1 , ■ ■ ■ ,PPi^ '■ 



3 The fractional stochastic convolution of the 
Poisson kernel 

In this section, we consider a bounded domain with C°° boundary, D C R fc , if 
k > 2, and D =]0, 1[ if k = 1. We consider the Green function defined in (2) 
and Gp(x, y) — (x A y) — xy, respectively. The purpose is to define the stochastic 
convolution j' D G k D (x,y)dB H (y) with respect to the fractional Brownian field with 
parameters Hi e [|,1[, i = l,...,k, introduced in the preceding section, and to 
study its sample paths. 

Throughout the section we shall make use of the following remark: Let k > 2 
and assume that for some norm || • || defined on a space of functions on R fc , we have 
sn Px£D\\G k (x,-)\\ < oo. Then, sup xeD \\E X (G k (B T , -)^J \\ < oo, and consequently, 

su PxeD II^d( x > *)ll < 00 • Indeed, since K x is a convex operator, denoting by P x the 
probability law of B T , we obtain 



E x (G k (B T , 



< E, 



G k (B T ,-) = f P x (dz) G k (z,-) <sup G k (z,-) 



We are interested in the integrability properties of G k D . To start with, let us state 
a result that for dimensions k > 3 is Lemma 2 in [20]. The extension to k — 1, 2 is 
trivial. 

Lemma 3.1 For any p E [1, mjkz^b there exists a positive constant JC\ depending 
on p and k, such that 

supllGkfoOHipp) </Cl (19) 

xeD 

For the values k = 1,2, (19) holds for any p e [l,oo[. Therefore by virtue of 
(18) we can choose p > sup ie | l5 ^ such that 



sup \\G D (x, -)|| i .... i < Csup \\G D (x, -)Hlpo(d) < oo- 

x&D L H i H k(D) x eD 



(20) 



Consider now the case k — 3. Property (19) holds for any p G [1,3[. Hence there 
exist po G [sup ig { lj k y jf, 3[ such that (20) holds. Then, accordingly with the results 
stated in the preceding section, for k — 1, 2, 3, T(G k D (x, •)) is a well-defined random 
variable in L 2 (Q) for any x £ D. 

A similar conclusion holds true for k = 4 under the additional assumption Hi G 
]|,1[, that is, excluding the possibility of having a standard Brownian motion in 
some of the components of B H . By similar arguments, the existence of T{G k D (x, •)) 
for k > 4 is ensured by the stronger hypothesis Hi G]^, 1[, for any i = 1, . . . , k. 
As we show in the next Lemma, one can relax this assumption by working with LP 
spaces with mixed norm. 

Throughout the section we suppose that D C [— R, R] k for some R > 0. We 
shall use the following inequality on Euclidean norms: For any /x > and $ > 0, 
z = 1, . . . , k, such that Yh=\ Pi — 1, 



i=i 



"<nw^ (2i) 



Lemma 3.2 Let A; > 4 and assume that X^Li Hi > k — 2. Then, there exists a 
positive constant /C 2 depending on H , D and k, such that 

su V \\G k D {x,.)\\ i ... i </C 2 . (22) 

Proof. By applying the inequality (21) with ft := 2 — k and $ := — , and the 
remark at the beginning of the section, we obtain 

/ fR / fR 2-k \ 

•>'.*•-*««, * 2 {L~iU x - y ^" m ) -*» 



< 2[] / ki-j/ii *i % 

i=i W-.R / 



The supremum on x G D of the last term is finite if and only if > — 1. By 

the definition of this condition is equivalent to X^ =1 Hi > k — 2. ■ 

In the sequel we will assume the hypothesis: 
(H) Hi G [|, 1[, i = 1, . . . , fc, and for dimensions A; > 4, we suppose in addition 
that Eti Hi > k - 2. 

Then, for any dimension k > 1, Lemmas 3.1 and 3.2 yield 

sup||^(x,.)|| i,., i </C, (23) 

with /C = max(/Ci, /C 2 ). 

This proves the existence of a stochastic process 

J = \j{x) := J (G* (x, •)) = j D G k D {x, y)B H (dy),x G , (24) 
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with values in L 2 (Q), satisfying 



sup E 

x&D 



f G k D (x,y)dB H (y) 2 = sup / (K*/ k) G k D (x,-))\y)dy 

D x&DJD v ' 



<C H ^ v \\G k D {x,-)f i .. . i <C H JC 2 <oo. (25) 



Next we prove that the stochastic field J has a.s. Holder continuous sample 
paths. 



Theorem 3.3 Under (H) it holds that 
G k D (x,-)-G k D (z,.) 



7r i''"'' 77 k (D) 



(26) 



/or any x, z E D, with 



A=l, 



A e 



for k — 1 , 



0,(2-HEti^)Al[, fork>2. 
Therefore, the Gaussian random field J defined in (24) satisfies 

E (1.70*0 -^)l 2 ) < C|rr-z| 2A , 



(27) 



and a.s., £/ie sample paths are Holder continuous of order 7 e]0, A[. Moreover, for 
any p E [1, oo[, J E LP (Q; L°°(D)). 

Proof. Fix x, z E D. For = 1 and D =]0, 1[, direct computations yield 

G k D (x,y)-G k D (z,y)\<C\x-z\, 

which implies (26). 
Let k > 2 and set 



2 



1T ^''"' 7r k (D) 



T« = ||G fe (x,-) -G k (z,-) 
fW = \\E x (G k (B T ,-))-E z (G k (B T , 

By the strong Markov property of Brownian motion, 

E z (G k (B T , •)) = E, (G k (B T - x + z,-)). 

Thus, similarly as in the remark at the beginning of the section, 



(28) 



f W < SU p G k (y, •) - G k (y -x + z, 

yeD 



L H i' 



Therefore, we can concentrate on proving that < C\x — z\ 2X for the values of A 
given in the statement. 
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By the very definition of G k , 



G k (x,y)-G k (z,y) 



< C 



\x-y\ 
z-y\ 



t^ k dt 



Fix A e]0, 1] and consider the change of variable t := (9\x — y\ x + (1 — 9)\z — y\ x ^j A , 
9 E [0,1]. Then, 



G k (x,y)-G h (z,y) 



< C 



x 



\x — y\ x — \z — y\ x 



f Q (9\x-y\ x + (l-9)\z-y\ x ) 



2-k-X 
X 



de 



< C\x- z\ x (\x- 



y\ 2 ~ k ~ x +\z - y\ 2 ~ k ~ x 



)■ 



(29) 



(30) 



where in the last upper bound we have used the inequality (a + b) p < dP + bP, valid 
for any a, b > and < p < 1. 

We can now apply Minkowski's inequality for L p spaces with mixed norm (see 
[3], page 302) to obtain 



T (k) < c i 



x — z\ 



|2A 



X 



x — 



12— fc— A 



777' "•73T 



+ 



L H i 



k(D) 



12— fc— A 



W "•73T 



*(£>) 



Suppose that A < 2 — A; + Y$=i Hi- Following the arguments of the proof of Lemma 
3.2, the inequality (21) applied to fi := A + k — 2 yields 

< C\x-z\ 2X . 

This ends the proof of (26) and, by the isometry property of the stochastic integral 
given in (25), we obtain (27). 

Since the process J is Gaussian, the statement about regularity of its sample 
paths follows from Kolmogorov's continuity criterion. More precisely, we get 

\J(x)-J(z)\ . 
sup - — := C(u) < +oo, a.s., 

x,zeD,x^z \x - zyi 

for any 7 g]0, A[ (see [24], Theorem 2.1, pg. 26). 

From here, we easily obtain sup :rgD [J7"(x)| < +00, a.s. This implies that J G 
L p (tt; L°°(D), by Theorem 3.2 of [2]. ■ 

Remark 3.1 Consider the following assumption which is stronger than (H): 

(H*) Hi e [|, l[ k , i = 1, . . • , k, and £f =1 Hi > k - 1 for any k>2, 

In this case, it can be proved that (26) holds true with A = 1 for any dimension 
k > 1. 
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4 Existence and uniqueness of solution to the 
fractional Poisson equation 

This section is devoted to establish the existence and uniqueness of solution to the 
equation (1). This result will be obtained by a pathwise argument; once it will be 
established, we will prove some probabilistic properties of the solution. We borrow 
the method of the proof from [20] (see also [4], [5] and [11]), which follows the 
classical monotonicity methods. We shall denote by C(D) the space of continuous 
functions on D and set S = {«; : w G C(D), w\qd = o}. 

For its further use we highlight some properties. The first one, denoted by (M) 
is a monotonicity property. The second one, named (P), has been proved in [4] 
(Lemma 2.4); it is a consequence of the solvability of the Dirichlet problem on D 
and Poincare's inequality (see [9] or [1]). They are formulated as follows: 

(M) / is a function of the form f = /1 + /2 with f\, f'2 : IR — > M, f\ non- decreasing 
and ji Lipschitz with Lipschitz constant L, if and only if for every u, v G M., 

(u — v)(f(u) — f(v)) > —L(u — v) 2 . (31) 

(P) There exists a constant a > such that for any tp G L 2 (D), 

J d (J D G k D (x,y)<p(y)dy^<p(x)dx < -a (J D G k D (x,y)<p(y)dy^ dx. (32) 

We begin with the existence and uniqueness result. 

Theorem 4.1 Assume (H) and suppose that g is a real function defined on D, 
1 ... 1 

g G L 1_H i' ,1 ~ H k( y D), f — fi + f 2 and satisfies (fl) and (/2) with a Lipschitz 
constant L < a. Then, there exists a unique stochastic process solution to (3) with 
sample paths in S, a.s. 

Proof. Consider the operator T : S — > S, defined by 

T(w)(x)=w(x)- I G k D (x,y)f(w(y))dy. 
Jd 

Holder's inequality for LP spaces with mixed norm and (26) yield 
(G k D (x,y) - G k D (z,yj) g(y)dy 



D 

< G k D (x,-)-G k D (z, 



4r,^ \\9\ 



L H i' ' H k(D) L 1 -"!' ' 1 - H k(D) 

<C\x-z\ x , (33) 

with A > given in Theorem 3.3. With this and the continuity of the stochastic 
convolution established in Theorem 3.3, we conclude 

b(x) := / G k D (x,y)g(y)dy+ [ G k D (x,y)dB H \y) G S, a.s. (34) 
Jd Jd 
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We next show that the operator equation Tw = b has a unique solution for any 
b G S, or equivalently that T is bijective. Uniqueness guarantees the measurability 
of the process {w(x), x G D}. 

The one-to-one property of T follows by applying (M) and (P) (see Theo- 
rem 2 in [20]). Indeed, let u, v be such that Tu = Tv. Then u(x) — v(x) = 
$ D G k D (x,y)(f{u{y)) — f(v(y)))dy. Multiplying both sides of this equation by 
f(u(x)) — f(v(x)), integrating over D and applying property (P), we obtain 

/ (u(x) -v(x))(f(u(x)) - f(v(x)))dx 

JD 

= J d (f(u(x)) - f(v(x))) G k D (x, y) - f{{y))) dy) dx 

< -a J d (/ d G k D (x, y) (f(u(y)) - f(v(y))) dyj ' dx 
= —a / (u(x) — v {x)) 2 dx. 
Because of (M), 

/ (u(x) - v(x)) (f(u(x)) - f(v(x))) dx>-L [ (u(x) - v(x)) 2 dx. 
Jd Jd 

Hence, (a — L) J D (u(x) — v(x)) 2 dx < 0. Since L < a and u, v G L 2 (D), this implies 
u(x) = v(x) for almost every x G D. 

We now prove that T is onto. In the next argument, we fix a sample path of the 
process B H on a set of probability one. 
Step 1: A solution for a regular problem. 

We follow the arguments of [20] Lemma 3 and the sequent discussion. Let b G S 
and b n G C^°(D), n > 1, be such that b n — > b in L 2 (D). Then, one can construct 
a sequence of functions solving Tu^ = b n and satisfying 



u 



(n) 



(x)= [ G k D (x,y)f(u^(y))dy + b n , for x G D, 

JD 



U 



(n) 



3D 



= 0. 



(35) 



By the properties (M), (P) and since w (n) G L 2 (D), one can prove that {u^ n >, n > 
1} is a Cauchy sequence in L 2 (D). Let u denote the limit. 
Step 2: u is the solution. 

We would like to pass to the limit (35) and obtain 



u 



\ x )= I G k D (x,y)f(u(y))dy + b(x), for x G D, u\ dD = 0. 

JD 



(36) 



By taking a subsequence of (still denoted by u^), we can assume that 
converges to u, a.e., as n — > oo. Then, since /i is continuous fi(u^) — > fi(u), a.e., 
as n — > oo. Consequently, Schwarz inequality with respect to the measure on D 
with density \G D (x,y)\dy implies 

G k D (-, y ) \fM n \y))-h{<y))]dy 2 



L 



D 

< c 



L2(D) 



dx / dyG k D (x,y) \fi(u {n) (y)) - fi(u(y))] dy. 



D 
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By bounded convergence, this last expression goes to zero as n — > oo. 
Similarly, 



/ G k D (-,y) [f2{u^\y))-h{u{y))]dy 

JD 1 J 

<C I dx I dy G k {x,y) \u^ n \y) - u 
Jd Jd 

< Csup \\G{x, OI|Li(£>)||« (n) - u\\h( D ), 



2 

L 2 (D) 



xeD 

which clearly tends to zero as n — > oo. 

Thus, u satisfies (36) and u E S. ■ 

The last part of the section is devoted to a further analysis of the solution u. In 
the next Lemma we set IC := sup xeD J D G D (x,y) dy. 

Lemma 4.2 Assume the assumptions of Theorem 4.1 and in addition, LJC < 1. 
Then, for any p G [1, oo[, 

IMlLP(n;L°°(D)) — 

Proof. The expression b given in (34) belongs to L P (Q; L°° (£))). Indeed, the function 
x — > J* D G D (x, y)g(y)dy, defined on .D, is continuous and deterministic and therefore 
belongs to the space LP (Q; L°° (D)) . As for J D G D (x,y)dB H (y), this property has 
been proved in Theorem 3.3. 

By definition, u(x) = f D G D (x, y)f(u(y))dy + b(x). Then, by the properties of / 
we have 

IMIlp(J2;L°°(D)) — & {l + L ll U llz,P(n;L°°(D))} j D Gd( X iV) + II &II Lp(SI;L°°(D)) • 

Taking into account the restriction on the constant L, this yields the result ■ 

We finally state a result on the regularity of the sample paths of the solution. 

Theorem 4.3 With the same hypotheses as in Theorem 4.1, for any x, z G D and 
for any p G [1, oo[, the solution u to (1) satisfies 



E(\u(x) - u(z)\) p < C\x-z\ pX , 

with A defined in Theorem 3.3 (see also Remark 3.1). 

Consequently, a.s. the sample paths are ^-Holder continuous with 7 g]0, A[. 

Proof. We write \\u(x) — u(z)\\ LP , n -, < X)f=i Iii x , z )i with 

= / (G k D (x,y)-G k D (z,y))f(u(y))dy 
Jd v 7 £P(n) 

/ (G^^^-G^z,?/))^^)^ 
/ (C&(x,y)-C&(*,y))d2?"M 

J D ^ ' lp (SI) 



(37) 



J 2 fe (x,z) = 
J 3 fc (x,z) = 
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By Holder's inequality and (26) we have 
lHx,z)< G k D (x,-)-G k D (z,-) 



i i 



x M+|/ 2 (0)|+L EH P 1 



k i D )/ 



<C [M+\f 2 (0)\+L [E\\u\\ p 1 



\x — z\ 



The factor Elk || p 1 , 

L 1 ~ H i ' ' 1 - H k(D) 

4.2. Consequently 



is finite. Indeed this follows from (18) and Lemma 



Ii(x,z) < C\x-z\ x 

In a similar but easier way, we obtain a similar bound for If (a;, z \ As for if (x, z), 
the bound is obtained by first applying the hypercontractivity inequality and then 
(27). 

This ends the proof of (37). The statement about the regularity of the sample 
paths follows from Kolmogorov's criterion. ■ 

5 Lattice approximations in L 2 -spatial norm 

This section is devoted to give finite difference approximation sequences for the 
SPDE (1) on the domain D =]0, l[ h , obtained by discretizing the Laplacian operator. 
To simplify the notation, we shall omit the superscript k when referring to G k Dl and 
denote by C any positive, finite constant not depending on n, although its value 
may change from one line to the next. 

The analysis of the rate of convergence is done by means of a general result given 
in the next theorem. 

Theorem 5.1 Assume the hypotheses of Lemma \.2. We consider a sequence of 
functions {g n }n>i defined on D such that 

lim \\g - g n \ i i = 0, 
along with another sequence of functions iGonf , defined on D x D satisfying 

l ' J n>l 



Id 



G D , n (x, .) - G D (x, Olj^,..,^ dx < Cn\ (38) 



ID 

for some 7 > 0. 

Let {u n (x),x G -D} n >i be a sequence of random fields satisfying 
u n (x)= G D , n (x,y)f(u n (y))dy+ G D , n (x,y)g n (y)dy + G D<n (x,y)dB H (y). 

JD JD JD 

(39) 
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We suppose that for some p G [1, oo[ and q G 



l-max(ifi,...,if fc ) 



oo 



SUp \\Un\\ L p(n;Li(D)) 
n>l 



< c. 



Then, 



\ u - u n \\ LP{n . L2{D)) < C in l + \\g-gr, 



(40) 



(41) 



where u is the solution of (3). 



Proof: We shall follow the scheme of the proof of Theorem 2.4 of [11] (see also 
Theorem 4 in [20]). By defining 

T{x) = JjG D {x,y) - G Din (x,y)]f(u n (y))dy 

+ / G D (x,y)[g(y) - g n (y)]dy + / [G D (x, y) - G Din (x, y)]g n {y)dy 
jd jd 

f [G D (x,y)-G D , n (x,y)]dB H (y), 
Jd 



+ 



we have u(x) - u n {x) = J D G D (x,y)[f(u(y)) - f(u n (y))]dy + T(x). 
Using (31) and (32), as in [20] we obtain 

(a - L)\\u - u n \\ 2 L2{D) 

<2a f (u(x) - u n {x))T{x)dx + f (f(u(x)) - f{u n {x))) T{x)dx. (42) 

JD JD 



Let q G [ i_ max( H lv „ ;gfc) , ooj and let q G [1,2] be its conjugate. By applying 
Holder's inequality and by virtue of the assumptions on /, the right-hand side (42) 
is bounded by 



(2a + L)\\u - u n \\ Lq ( D) \\T\\ L ^ D) + 2M||T|| £S 



?(£»)■ 



Lemma 4.2 yield 



sup \\u u n \\ LP , n . Lq , D ^ < C, 

n>l 



(43) 
(44) 



for the values of p and q such that (40) holds true. 

From (42)-(44) and applying Schwarz's inequality we obtain 



\u — u r 



<C \\T 



\Lv(Q.;L i i(D)) 



\LP(Q;L 2 (D)) 

Holder's inequality for LP spaces with mixed norm yields 

|T(x)|< \\G D (x,.)- G D , n (x, -)|| fn JI/WII tT ^- 

+ \\G D (x, -)|| JL-,...,j- t \\g-9n\\ rj — Kj—,--- ,j 
L H i H k(D) L 1 ~ H i '! 



(45) 



(D) 



+ \\G D (x, •) - G D:Tl (x 

+ 



~ H k (D) 
9n\\ T - 



D 



L H i' ' H k (D) L 1 ~ H i ' ' 1 - H k(D) 

[G D (x,y)-G D , n (x,y)]dB H (y) . 
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By the assumptions on the function /, we have 

\\f^)\\ L ^...^ D) <M+\hm+m 

Thus, by (18) and (40) 

sup 



n>l LP(n;L 

Consequently, since q < 2, 



\T\\LP(n-,L^(D)) < C\\T\\ LP (n. L 2( D ^ 



< c. 



*(£»)) 



< C 2 l + 



+ 



f D [G D (-,y)-G D , n (-,y)]dB H (y) 



LP(Q;L 2 (D)). 



(46) 



where we have applied (23) and 7 is given in (38). 

Let us now give an upper bound for the last term in (46). Assume first p G [1,2]. 
By applying Holder's inequality with p — | > 1 to the expectation operator and 
then Fubini's theorem, we obtain 

/ [G D (;y)-G D , n (;y)]dB H (y) 
Jd 

E I [G D (x,y)-G D , n (x,y)]dB H (y) 
< [E I [G D (x,y)-G D , n (x,y)]dB H (y) 



LP(Q;L 2 (D)) 



D 



L2(D) 
2 

L 2 (D) 



< 



< 



D 



D 



E 



D 



[G D (x,y)-G Dtn (x,y)]dB H (y) 



dx 



G D (x, •) - G D , n (x, •) 



1 .... 1 dx 

L H i H k(D) 



< Cn-1, 



where in the last inequality we have used (38). 

Next, we consider the case p > 2. We apply first Minkowski's inequality with 
respect to the probability measure and the Lebesgue measure, then the hypercon- 
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tractivity inequality (see for instance [18]), to obtain 



D 



[G D (;y)-G D>n (;y)}dB H (y) 



LP(H;L2(D)) 



E 



D 



D 



[G D (x,y)-G D , n (x,y)]dB H (y) 



dx 



f f [G D (x,y)-G D , n (x,y)]dB H (y) 
Jd Jd 

- ( L II ( L [Gd{x > y) ~ y^ dBH ^ 

= ( J d (e| jf [G d (x, y) - G D , n (x, i/)]dfl*(y) 



zJ(f2) 



<(p-l) 5 



E 



dx 



[G D (x,y)-G D , n (x,y)]d5"(y) 



2\ 



jdx 



D 



G D (x, •) - G DjTl (x, ■) 



L^'"''^(D) dX 



<C(p) 
< C(p)n~2. 

Plugging these estimates in (46), we finish the proof of the theorem. ■ 

Let I k and J* be the sets of indices {1, 2, . . .} h and {1,2,..., n — l} fc , respectively. 
On the space X = {u : w = {ui} ieI k} = R^ 1 ^ endowed with the Hilbert-Schmidt 
norm, we define the second order difference operator A : X — > X, 

k 

(Au)i = J2 n2 ( U i~ej ~ 2u i + U i+ej ) , i £ 
J'=l 

where {ej}^ =1 is the canonical basis of M. k . 

Consider the orthogonal complete system in L 2 (D) provided by the functions 

vp{x) = sin(/?i7rxi) • • • sin(/3 fc 7T£ fc ), (3 (z I k , k > 1. 

The set of vectors |^|^ ^"/3> e ^n}> (^s)i = ^(n)' * e ^ ^ s an orthonormal 
system in X of eigenvectors of A, with eigenvalues \p = —^(fifc^ + h /^C/jJ, 



where q = ( ) sin^ ( tt ) • Notice that -% < Q < 1 for every 1 < / < n - 1. 



2ri 



2 / Ztt 



2;( 



We will consider approximation schemes of (3) based on the grid of D given by 
Q 





(-•■ 


3k\ 


In 


V n 


■' n) 



For any point ^ G we set Dj = [^-, 
define K n (x) = 1 . 



x • • • x 



rife 2fc±I[ and for each x G D j we 
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We begin by giving a first type of discrete approximations of u on points of Q, 
denoted by u n . If ^ G Q fl 3D, we set u n (£) = (boundary conditions), while for ^ 
with j G J*, we define « n (^) to be the solution of the system 

^n = /K)+4n + ^B", (47) 

where B H is the vector {B H (A) = / R * t Di {y)dB h r (y), z G /*} and (a n ) n >i a 
sequence of step functions, g„(x) = n > 1. Then, for any x E D we 

define u n (x) = «„(«„(#)). From [11] we know that {ii n (x),x G -D} satisfies the 
evolution equation 

Un(x)= / G Dtn (x,y)f(u n (y))dy+ G Dtn (x,y)g n (y) + G Dtn (x,y)dB H (y), (48) 



where 

2 fc 

Gr>,„(x,j/)= J! -r-^(«n(^))^(«n(y))- (49) 

pel* X P 



In dimension k = 1,2,3, we shall consider {u G D}, n > 1, as sequence 

of approximations of the process {u(x),x G D}. We notice that in this case the 
kernel Gd,u(x, •) is related with the truncation of the Fourier expansion of Gd{x, .)• 
For k > A we shall follow the more sophisticated approach of [20], which considers 
a smoothed version of Gd(x, .). We remark that for such dimensions G\d(x, .) is not 
square integrable. 

For low dimensions, we have the following. 

Theorem 5.2 Assume the assumptions of Lemma Let k G {1,2,3} and 

{u n (x),x G D}, n > 1, be defined in (48), (49). Then, for any p G [l,oo[ ; 

\\u - u n \\ LP{n . L 2 {D)) < C I n~ v + \\g - g n \\ 2 _^_ ... ], (50) 

V L 1 -"!' ' 1 ~ H k (D) / 

with v G]0, |], ^ g]0, \[, v g]0, |], /or k = 1, k = 2 and = 3, respectively. 

Proof: It follows from Theorem 5.1. Indeed, first we give a slight improvement of 
Lemma 3.4 in [11] which yields the validity of condition (38) for Gd, u '■= Grj, n with 
7 := Av and the values of v of the statement. 

In fact, for the expression termed A in [11], page 223, we have 



roo 

A<C r k - l - 4 dr = Cn k ~\ 

J n 



while for k = 3, we can proceed with the term called B as follows. Let p G]l,2[, 
then 

1 C 1 



. ,a\ 2 n 2 nP ^-t k \a\ A ~P 



< Cn-\ 
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The process u n := u n satisfies (40). This can be easily checked by applying 
Lemma 3.3 in [11] and similar arguments as those in the proof of Lemma 5.7 below 
(we leave the details to the reader). ■ 

We next deal with higher dimensions. The Fourier analysis techniques we 
shall use in the proofs require the identification of functions / : [— 1,1 [ fc — > K 
with functions F : T k — > H. defined on the k-th dimensional torus through the 
exponential mapping F(e l7TX ) := F(e tnxi , . . . ,e t7TXk ), which carries Lebesgue mea- 
sure into the Haar measure, that is jj-i^pt f(x)dx = j T k F(e inx )dx. When deal- 
ing with the function y — > Go(x,y), we will consider its odd extension, that 
is, for any x,y G]0, l[ fe , y = (y u ...,y k ), we define G D (x, (y u . . . , -y u . . . , y k )) = 
—Gd(x, (y±, . . . , yi, . . . , yk))- We still note Gd(x, •) the extension. Let <Gs x D (e my ) = 
Go(x,y) be its identification with a function defined on T h . Observe that Gfj sat- 
isfies 

||G£(e i7r -)l|L^-^(T*) = \\G D (x,-)\\ LPl ,..., PH] ^ m = 2 E ^^\\G D (x,-)\\ LP1 ,.., Pk{D) , 

for any pi, . . . ,p k such that the last norm is finite. 

Let ip(x) G C~(] - 1,1[) be an even function, < ip < 1, f\ ip = 1. Set 
^(x) = IljLi i'ixj)- Clearly, fy(x) G C£°(] - 1, l[ k ) and it is an even function in each 
variable Xj. Define 

$( e *™) = fl 0(e^) := f[ i>( Xj ) = *(x). 

3=1 3=1 

The functions <& e (e mx ) = pE-^(f) := ^(z), £ > 0, provide an approximation of the 
identity in T k . 

We shall denote by \1/ the Fourier transform of VP, which is a rapidly decreas- 
ing function, therefore for any 9 G [0,oo[ there is a constant C{6) such that 

Bu Pc iefl*(fli<c(0)- 

Let us now introduce a second kind of approximations of u. For this we start by 
writing A = U t AU, with U the (n — l) k matrix whose rows are the vectors U^, (here 
(3j, j = 1, ■ • • , (n — l) k is the lexicographic enumeration of I k ) and A the square 
diagonal matrix with entries Ajj = 

The smoothed version of A is defined as follows. Fix e > and define A e as the 
square diagonal matrix in dimension (n — l) k with diagonal entries 

A £ = A/3j 
ft 

In connection with A £ we define a sequence (u* , n > 1) of functions in the following 
way. If £ G QndD, set m^(^) = (boundary conditions). For with j G define 
M n(n) to be the solution of the system 

(E/*A e Z7K = /«) + + n k B H . (51) 
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Finally, for any x G D we define u £ n (x) = u £ n (n n (x)). 

We shall prove later that an appropriate sequence u n := u 6 ^ of such functions 
converges to the solution of (1) in the space LP (Q; L 2 (D)) , for any p > 1, with a 
rate of convergence which depends on the dimension k, on the driving noise and on 
the rate of convergence of g n to g. 

The following result is proved with the same arguments as in Proposition 1 of 
[20]. 

Proposition 5.3 With the same hypotheses as in Theorem 4-1 and assuming that 
the Lipschitz constant satisfies L < 4k, Equation (51) possesses a unique solution. 
Moreover, this solution satisfies the mild equation 



(x) = J D G^ n (x,y)f(u £ n (y))dy + J D G t Di Jx,y)g n {y)d V + G £ Dn (x, y)dB H (y), 

(52) 

where 

^ i>{ej3)2 k 

G e D , n (x,y) = J2 — \ vp(K n (x))vp(K n (y)). (53) 

f>eik Af3 

Both (49) and (53) correspond to discretized Fourier series expansions; in (53), 
the Fourier coefficients are smoothed by the factor ^(ef3). 

Our next aim is to apply Theorem 5.1 to u n := u^ n ^ defined in (52), for values 
of e that depend on n, and dimensions k > 4. The next statements provide the 
ingredients for checking condition (38) for Go,n '■= G"d^- 

Set G £ D (x,y) = G^ £ (e^), where G x ^ £ (e iwy ) = J Jk G x D (e i7T{y ' u) )^e(e iwu )du. The 
function G e D (x, y) is a smoothing of Gd(x, •). 



Lemma 5.4 For any e > 0, we have 



-i>{e(3)2 



G £ D (x,y)= n 2^2 M x )Mv)> ( 54 ) 



in L (D x D) and a.e. In addition, 



ll^(x,OII| 2(D) = ^E^^^) ) (55) 
and the series converges uniformly in x e D and e e]0, Eq\. 

Proof. The first part of the assertion is Lemma 8 of [20] . For the sake of completeness 
and further use, we give some details of its proof. Fix p e 1, • Young's 
inequality for convolution ([1], page 34, Corollary 2.25) yields 
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\GUx r )\\ LHD) = 2-^\\G^(en\\L H ^) 

<2 2 ||Gf ) ||£ P (Tfc)||^e||L'-(T fc ) 

= C\\G D (x, •)||LP(D)||^ r e||L'-a-l,l[fc), 



for \ = i + i - 1. 

2 r p 

For any r > 1, sup £ ||^ r e ||L'-(]-i,i[fc) < C- Thus, using (19) we obtain 

sup sup \\G £ D (x,-)\\ L 2 {D) < C. 

xeD £ £]o, £0 ] 



(56) 



Hence G £ D £ L 2 (D x D) and the formula (54) follows from the computation of the 
Fourier coefficients carried out in [20], Lemma 8. 

The orthogonal complete system (vp,/3 G I k ) satisfies ||v/3||l 2 (d) = 2~ fc / 2 . Thus 
(55) follows easily from (54). Finally, (56) implies the uniform convergence of the 
series in (55). ■ 

The next result provides an estimate of the discrepancy between Go and G £ D . 
Lemma 5.5 Assume (H). Then for every e > 0, 

< Ce\ 



sup \\G D (x, •) - G e D (x, 

x&D 



777 Hi, 



with 



A=l, 



k(D) 



for k — 1 , 



(57) 



A e 



0,(2-HE-=i^)ai[, fork>2. 



Proof. Since Gd{x, •), G £ D (x, •) are odd in the ^-variables, we have 

\\G D (x, •) - G £ D (x, -)IU,.., im , = C\\G D {x, •) - G £ D (x, 0L 1 ^ 



= C\\G x D (e iw ') -G^e* 



with C = 2 Si=i H \ Using that J T * <E> e (e"™)cfot = 1, we can write 



L H i H fe (T fc ) 



/ (G^(e m ') - G x D (e iw( '~ u) ))$s(e i7TU )du 

J7 k 



l i 
L H i' "' H k (T fc ) 



< 



G*( e iT(-«))|| ^ ^ $ e (e l ™)d« 



= 2£*=i J * / ||G D (a;,.)-Gi ? (a;,.- 

<cf 



u 



i i 



L i 



$ e (e l ™)Ai 



M A $ £ (e"™)(2w 



C 



/ 

vi-i,if 



e k 



x 



22 



where we have applied Minkowski's inequality with respect to the finite measure on 
T k defined by § £ (u)du and eventually (26). From this we infer (57) by observing 
that the support of the function ^ is included in ] — 1, l[ k . 

■ 

As an additional auxiliary result, we need a priori estimates for the solution of 
(52). An ingredient for this is provided by the following Lemma. 

Lemma 5.6 Fix e > 0. The smoothed, discretized Green function defined in (53) 
satisfies 

sup sup \\G e Dn (x, -)IU 2 (D) < +oo. (58) 

n>l x£D,ee]0,E ] 

Proof. The system {vp{K n {y))} is orthogonal in Rt™ -1 )*, thus in L 2 {D) as well. 
Hence, using the lower bound jA^] > A\f3\ 2 we have 

< Mi2 ^ 2 (^)2 fc 2 , . 2* &{e(3) 2 

From the last inequality, along with (55), we obtain 

sup sup \\G £ D (x,-)\\ 2 L 2 {D) <C sup { E ^\mP v U x ) 

n>l x£D,e£]0,e ] xeD,ee]0,s o ] ^ Q j k \fJ\ 



I3ei k 

<C sup \\G £ D (x,-)\\ 2 L2(D) < +oo. 

xeD,ee]0,e o ] 



Lemma 9 of [20] gives a more particular statement than the previous Lemma 
5.6. We have found an incorrect argument in the proof of the former that can be 
fixed using the proof of the later. 

Let q G [1,2]. By (18), Holder's inequality and Lemma 5.6 we have 



sup sup (\\Gh tn (x,.)\\ V||G^ n (av)|k(D) J < /C, (59) 

n>l xeD,ee]Q,e ]\ L H i H k (D) J 

for some positive, finite constant K. 

We can now prove an a priori estimate for the solution of (52). 

Lemma 5.7 Assume the same assumptions as in Proposition 5.3 with the Lipschitz 
constant satisfying the restriction L < min(4/c, JC^ 1 ), where JC is given in (59). 
Then, for any p G [1, oo[ and q G [2, oo[, 



sup sup (\\u £ n \\ LP( n. Lq(D)) ) < C. 

n>l ee]0,e ] V 
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Proof. Fix q G [2, oof and denote by q g]1, 2] its conjugate. Holder's inequality and 
the properties on / imply, for any x G D, 

|<(x)|<sup sup \\G%(x r )\\ L ^ D) (M+\f 2 (0)\ + L\\u' n \\ Lq(D) ) 

n>l x£D,s£]0,en] \ / 



x£D,e£]0,e 

sup sup \\G £ Dn (x 

n>l x€D,e€]0,£ ] 



i .... i , , SUp n 



L 1T i' '' 1T k{D) n >l 



i-H. < — <T=lT 



HD) 



+ 



D 



G £ D Jx,y)dB H (y) 



(60) 



Since u £ n is a step function, its L 9 -norm is finite. Moreover, arguing in a similar 
manner as we did in Theorem 5.1 to get an upper bound on the last term of (46), 
and applying (59) we obtain 



D 



G%J;y)dB H (y) 



LP(n-,Li(D) 



<C(p,q,H). 



Taking the LPiVt; L <? (D)-norm in (60) yields the conclusion. 

The next lemma gives an estimate of the discrepancy between G £ D and G £ D n . 



Lemma 5.8 Assume k > 4. Let 5 G]0,2[ ; fi G 0, ^ and set e(n) = n M . There 
exists a positive, finite constant C(5,k), depending on S and k but not on n, such 
that 

\\Gf ] - G^\\l H d, D) < C(5,k)n- l >. (61) 

Proof. We follow the proofs of Lemma 3.4 in [11] and Lemma 10 in [20]. By the 
definitions of the kernels G £ D n and G £ D given in (53) and (54) respectively, we have 



Gd^ ~ G £ D™n\\h{DxD) < C^2Ai(x,y), 



with 



A 1 
A 3 



= 1 E — 

JDxD zri. 



2 k ^(e(3) 
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vp{x)vp{y) 



dxdy, 



* 2 W h 



2 k i>{e(3)vp{x)v f3 {y) 



dxdy, 



2 k i>(e(3) 



, A/3 



[vp(x) -vp{K n {x))]v p {y) 



dxdy, 



DxD 



E 



2 k i>(e(3) 



vp{n n {x)) [vp{y) - v p {n n {y))\ 



dxdy. 
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In the sequel, we shall write e instead of e(n) for simplicity, and we fix 5 > 0. 
Remember that (vp,(3 G I k ) is a family of orthogonal functions in L 2 (D) with 
Ik/3|U 2 (D) = 2~2. Let 9 > ^=^. Since & is a rapidly decreasing function, we have 



m= E 



&(e[3) ^ C{9) 



< 



7r" 



t29 



E 



i 



14+20 



< C(^) £ - 2e n^ 2e+fc = C(fl)n 2 ^~ 4 - 2e+A; . 

Fix := <5 ^ fc 2 ~ 4 in the last expression. We obtain A 1 < C(9)n~ s . 

For the analysis of the term A 2 we apply the estimate 
for any (3 G I k . Taking 9 < we obtain 



A/9 



< 



C 



, valid 



A, 



< 



E 

pel* 

m 

£ 28 n 2 



E 

Pel* 



1_ 2 



* 2 (^) < ^ v i!M 



2+29 



Pel* 

< C(6)n 2 ^ +k - 2e . 



Consider 9 : = f^g. The last estimates yield A 2 < C(9)n~ 5 . 

For the study of the remaining terms, we use that for any j3 G I k , \vp(x)—vp(z)\ < 
C\j3\ \x — z\, and \Xp\ > 4\/3\ 2 . This ensures 



. . a \ C ^ y 2 (eP)\(3\ 2 C{9) ^ 
max(A 3 ,A 4 ) < - E ^ 2 < 4f^ 2 E 
n ^e/* £ n Pel* 
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Hence, as for A 2 , we obtain max(A 3 , At) < C(9)n~ 5 . 

The proof of the Lemma is now complete. ■ 
As a consequence of Lemma 5.5, (18) and Lemma 5.8 we obtain the following 

result. 

Corollary 5.9 With the same assumptions as in Lemmas 5.5 and 5.8, there exists 
a positive constant C not depending on n, such that 



j \\G D (x,.)-G^(x, 



dx < Cn 



-7 



(62) 



with 7 = (2/xA) A 5. 



Assume (H*). Then in the preceding Corollary, A = 1 and 7 = 2/j, A 5. The 
biggest value of 7 occurs when 2/i = 5. Since ji < and the equation 2<y lZ 2 — $ 



has the solution 8 — |, we conclude that 7 G 



Suppose next that k — 1 > X^ fc =1 Hi > k — 2. In this case 

(4-2£) (2-fc + Eii^) 



2^A G 



0, 



k-2 
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As before, the biggest upper bound of 7 is obtained by solving the equation 



(4-2tf)(2-fc+X)jLi«i 



7G 



0, 



fc-2 
4( 2-fc+ 



5, which leads to the value 5 



4 ( 2 -*+Eti g <) 
2 - fe + 2 Eti^ 



Therefore 



Consequently, under the assumption (H), we obtain 



7 G 



0, 



4 (2 - A; + (£* =1 ifr) A(fc-r 
2-fc + 2((Eti^i)A(fc-l)) 



We have now the ingredients to establish the convergence of the discretized 
scheme defined in (52) to the solution of (1) when the parameters n and e are 
related by the same constraints as in Lemma 5.8. The constant /C in the next 
statement is given in (59). 

Theorem 5.10 We assume the hypotheses of Lemmas and in addition, L < 
min(4A;, JC^ 1 ). Let k > 4 and e{n) be as in Lemma 5.8. Then for any p G [l,oo[, 



u 



u 



e(n)\ 



with V G 



0, 



LP(n;L 2 (D)) ^ 
2 -*+fe=i H *W-l) 



< C[n- v + \\g-g n \ 



L 1 ~ H i ' 



(63) 



2-k+2 Eti^ A(fe-l) 



Proof. It follows from Theorem 5.1, Corollary 5.9 and the discussion that precedes 
the statement. ■ 



Let po := maxj=i v .. i )t (i^7r) anc ^ assume that 

||<7 — 9ti\\lpo(d) < C'n \ 



(64) 



with a constant C independent of n. Then the right hand side of (50) and (63) can 
be replaced by Cn~ u . We end this article by giving some examples of sequences 
{g n )n>i fo r which (64) holds. 

Let g n (x) = g{n n {x)). Consider firstly the case k — 1, and assume that g is a 
continuously differentiable function defined in ] — 1, 1[. Clearly, for any x g]0, 1[, 
\g(x) — g(n n (x))\ < swp^ <1 \g'(x)\ n~ l . Consequently, (64) holds. Consider next 
the case k > 2. Suppose that the function g is differentiable and Vg belongs to the 
Sobolev space W m ' Po (D) consisting of weakly differentiable functions up to the order 
m, with weak derivatives in L P0 (D). Assume m > ^. By the Sobolev embedding 
theorem (see [1], page 85, Theorem 4.12), Vg is a continuous function on D, and 
we also get (64). (See [11] for similar remarks for p — 2). 

Let g G W 1 ' Po (D). Define 

g n (x) :=n k Y, [J 9(y)dy) t Dj (x), x G D, 
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where I k — {0, . . . , n — l} k . Applying Equation (7.45) in page 157 of [9], we have 



/ g(x) -n k g{y)dy dx < (u k n^ 

JDj JDi v 



n 



J J D . 



\(Vg)(y)rdy 



<C(k, Po )n-™ (Vg)l Dj 



PO 



Since 



we obtain (64). 



\n - a \\ po 

\iJ yn\\LP0(D) 



= E 



D, 



g(x) -n k g(y)dy 

JDi 



PO 



dx, 



Remark 5.1 Assume (64)- By applying Borel-Cantelli's lemma, as in [11], we can 
prove the following statements: 

1. With the same assumptions as in Theorem 5.2, there exists an a.s. finite 
random variable £ such that 

\\u - u n \\ L 2 {D) < ^n~ u , 

a.s., with v e]0, \[ for k = 1,2, and v e]0, |[ for k = 3. 

2. Assume the hypotheses of Theorem 5.10 and let v be as in this theorem. Then, 
there exists an a.s. finite random variable £ such that a.s. 



\u-^\\ L ^ D) <in\ 
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